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1 Introduction 

The Grcgory-Laflamme instability [1, 2] has been explored principally via linearized 
perturbation theory and thermodynamic scaling arguments. In its simplest incarna- 
tion, it is a tachyonic mode in perturbations around an uncharged black brane which 
makes the horizon bulge out at some points and squeeze in at others. It is a purely 
classical effect in Minkowskian signature, and we may think of it crudely as a black 
brane redistributing its mass in order to gain entropy. Until recently it was commonly 
supposed that this evolution would proceed until the horizon bifurcated into an array 
of black holes (and then these black holes might run into one another, clustering mass 
into ever fewer, ever larger black holes). 

This naive view has recently been called into question because of a result of Horowitz 
and Maeda [3]. These authors demonstrated rigorously that a black hole horizon cannot 
pinch off in finite affine time, and they give less airtight arguments that such a pinch-off 
also cannot happen in infinite affine time. Thus it becomes plausible that the endpoint 
of the evolution is a stable, non-uniform black brane. ^ In the case of a black string, 

^Something not too distant from this picture was suggested in section 3 of [4], where it was also 
suggested that non-uniforniity would develop through a first order transition. The reasoning there 
was rather different, and much less rigorous. 
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the non-uniform solution might be expected to look like an unbroken string of identical 
beads. 

Briefly, the aim of this paper is to look for such a non-uniform black string using per- 
turbation theory around the uniform solution, and to discuss possible phase transitions 
between different types of solutions. 

The organization of the paper is as follows: in section 2, wc discuss order parameters. 
In section 3, we perform an explicit perturbation analysis around the black string solu- 
tion in five dimensions. Non-uniform solutions are constructed perturbatively around 
the uniform string with rj — rjc- Unfortunately, even this involves some numerical work, 
which will be explained in some detail in sections 3.2-3.6. In section 4 we make some 
general remarks about how the Landau theory of phase transitions might bear on non- 
uniform branes of various dimensionalities, and also comment on possible instabilities 
which could arise after a brane has become non-uniform. We summarize our key results 
on the five-dimensional uncharged black string in section 5. 

2 Order parameters 

The Gregory-Laflamme instability is an infrared effect: when it is present, a// accessible 
wavelengths below a certain threshold are unstable. Life is simpler if one imposes an 
infrared regulator: for example, one may compactify the directions parallel to the 
black brane on some torus. Then the instability shows up only if the Schwarzschild 
radius is sufficiently small compared with the size of the torus (the zero-mode cannot 
be excited because of energy conservation). The commonly supposed endpoint of the 
evolution (before [3]) was a single black hole, localized in the torus directions as well as 
the orthogonal directions. By a simple scaling argument, such a black hole has more 
entropy for a given mass than the uniform black brane starting point, provided the 
size of the torus is sufficiently large compared to the Schwarzschild radius of the black 
brane. 

As long as we are discussing solutions of the classical empty space Einstein equations, 
Rfiu = 0, there is no intrinsic scale in the problem. This may be seen from the fact that 
a rigid rcscaling of the entire metric takes solutions into solutions. The uniform black 
brane determines a length scale, namely its Schwarzschild radius, -Rschwarzschiid- After 
compactifying the uniform black brane on a torus, suppose the longest wavelength 
perturbation accessible has wave number /Cmin- Then we may define the figure of merit 

— -fi-SchwaTzschild^min ■ (1) 

For a non-uniform solution, we may still define r) as the value that would pertain if the 
non-uniform brane were replaced by a uniform one of the same mass. 
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If the black brane does become non-uniform, there is another dimensionless number, 
A, which describes how non-uniform it is. Supposing that the rotational symmetry in 
the directions transverse to the brane is preserved, we can define unambiguously the 
area of a slice of the horizon over any point along the brane, and hence an effective 
Schwarzschild radius at any point. Then we can define 



Thus A = for the uniform black brane; A = oo for a black hole localized on the torus, 
and A is finite for a non-uniform black brane with an unbroken horizon. 

A final meaningful number for a non-uniform horizon compactified on a torus (or 
other manifold for that matter) is the number n of local maxima of the Schwarzschild 
radius — again on the assumption that transverse rotational symmetry is preserved, so 
that the Schwarzschild radius has a good local definition everywhere on the brane. 

Roughly speaking, the problem gets harder as r] decreases. For rj large, the compact- 
ification scale is much smaller than another scale in the problem, so we may expect the 
uniform black brane to be stable and to win out over any other solution. Below some 
critical r)c, the uniform black brane becomes unstable to small perturbations: this is 
the Gregory-Lafiamme instability. What happens for smaller rj is not known. Focusing 
first on the black string, one might reasonably expect that for some range of rj below 
r/c, there is a stable solution with finite A and n = 1. Then the question becomes, 
docs non-zero A develop smoothly or suddenly as one crosses ?7c? That is, is the phase 
transition from uniform to non-uniform black strings continuous or first order?^ The 
results of section 3 will suggest that it is first order. 

3 Numerics on the black string 
3.1 Generalities 

At 1] = ?7c, the perturbative Gregory-Lafiamme instability becomes a zero-mode of 
the uniform black brane solution, so we may expect that there is a new branch of 
the solution space heading out to nonzero A, but emanating from A = and rj = rjc- 
This new branch has one maximum for the Schwarzschild radius {n — 1) because the 
zero mode changes it!schwarzschiid, multiplicatively, by an amount 1 -|- Xcoskx where k 
is the fundamental wave-number of the black string.^ A continuous moduli space may 

^It makes sense to speak of phase transitions even though the system is of finite extent because it 
is much larger than Planck scale. 

■^We defined the 1/2 in (2) so that it would coincide with the usage of A in this section, to lowest 
order in small A. We will not be concerned with a more precise translation of definitions. 




(2) 



3 



seem at odds with first order behavior; see however figure 1(a), equation (33), and the 
surrounding discussion for a preview of how the two are consistent. 

We want static five-dimcnsional black string solutions to Einstein's equations in the 
absence of matter: R^i, = 0. A general ansatz for such solutions is 



-e^^ ( 1 



2M 



df + e 



2B 



dr^ 



1 - 



2M 



+ dz' 



(3) 



where A, B, and C are arbitrary functions of r and z. Here we have used difTeomor- 
phism freedom to constrain the ratio of grr and g^z, to get rid of any g^z term, and 
to ensure that the black string horizon is at r = 2M. The ansatz (3) is almost in 
conformal gauge: it would be exactly in conformal gauge if we changed coordinates to 
f such that df^ — dr'^/ (1 — 2M/r).^ Clearly, if A — B — C — 0, we have the uniform 
black string. 

It's possible to change coordinates and rescale the entire metric so that it becomes 



ds' 



-e^^ 1 



df + e 



23 



l-i 

y 



+ dx^ 



+ e 



(4) 



The Schwarzschild radius of the uniform solution is 1. The periodicity of the x di- 
rection is 27r/?7, so vj coincides with the minimal wave- number k which we can excite. 
Looking for solutions of the form (4) to R^v — amounts to an elliptic boundary value 
problem: boundary conditions are imposed at the horizon to ensure that it is regular, 
and other boundary conditions are imposed at infinity to ensure that the solution is 
asymptotically flat. (Actually, we will see that regularity plus asymptotic flatness are 
not quite sufficient to fix all integration constants. This subtlety will flgure promi- 
nently in sections 3.4 and 3.5, but may otherwise be glossed over.) Elliptic boundary 
value problems for ordinary differential are fairly straightforward to deal with: the 
basic strategy is some form of "shooting" algorithm, where one artificially imposes full 
Cauchy data at one end of the integration region, then integrates the ODE's to the 
other end, and then checks to see whether the boundary conditions are satisfied there. 
To understand why elliptic boundary problems are harder for partial differential equa- 
tions, imagine trying the same strategy, for example imposing Cauchy data at y = 1 
and then integrating out to infinity. One has to discretize the x-interval to perform this 
integration. Intuitively, such a discretization involves about as many Fourier modes as 
there are pixels. The higher Fourier modes have higher mass. Each of these modes has 
a solution which grows exponentially at infinity rather than vanishing. Round-off error 
will tend to source these growing solutions, with the result that numerical integration 

^We are able to constrain the horizon to He at a fixed value of r after fixing the conformal gauge 
because this gauge is preserved by conformal transformations. An appropriate conformal map will 
indeed send the region outside the horizon to the half-plane f > 0, or r > 2M. 
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outward in y quickly diverges, for any initial data at the horizon. It appears in fact 
that there is no standard algorithm for solving elliptic boundary value problems for 
non-linear PDE's. A custom algorithm is usually constructed, based commonly on 
some global relaxation method, to deal with a specific problem. 

Instead of pursuing the intensive numerical route, let us make an expansion around 
the uniform string: 

A = A^ao(y) + Xai{y) coskx + X^a2{y) cos2kx + O(A^) 

B = A^6o(y) + A6i(y) cos kx + X%2{y) cos 2kx + O(A^) (5) 

C = A^co(y) + Aci(y) coskx + A^C2(y) cos2A;x + O(A^) . 

We will work at small A, but the coefficient functions are finite. The expansion (5) is 
appropriate for studying perturbations at the wavelength which is marginally stable. 
The reason is that at lowest order, 0(A), a perturbation at this wavelength is a zero 
mode of the uniform black branc. At O(A^), back-reaction leads to non-trivial Xq and 
X2, where X = a, b, or c. If we make a Kaluza-Klein reduction in the x direction, then 
the Xi for i 7^ are massive modes, while the Xq are massless modes. Thus the Xj 
for i ^ fall off exponentially for large y, while the Xq fall off as inverse powers of y. 
So to read off the mass density of the black string, all we need is the x-independent 
modes. We will elaborate on this point later. 

Having made the expansion (5), we are left with a finite set of ODE's at each order in 
A. Boundary conditions still must be imposed both at the horizon and at infinity. Let 
us now enter into the details of the calculation.^ The reader uninterested in technical 
points may wish to skip directly to section 3.7. 



3.2 Solution to 0(A) 

The basic ODE's can be derived by plugging (5) into the Einstein equations, Rf^^ — 0. 
At 0(A), one can solve algebraically for 61 (y), with the result 

a^ + 2 {-1+y) {2ci + ya[ + 2yc[) 
3^y ' 

where primes denote d/dy. This equation incorporates the zero energy constraint. One 
may also show that the 0(A) equations imply 

2 

b[ = _^^^y {k^y^ai + 2b,-2ci + 2k''y^ci + 2ai-2ya, + c,-2yci) , (7) 



^Readers interested in a fuller description can examine the Mathematica notebook [5]. 
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which is a relation which we will find helpful in section 3.4. Eliminating bi from the 
other differential equations at 0(A) leads to 



Evidently, the horizon is a regular singular point of these differential equations, which 
are linear because at 0(A) we're dealing precisely with linearized perturbations around 
the uniform black string. We can exploit the linearity to get rid of one arbitrary 
constant: let us say Ci(l) = 1. Examining the region near y = 1, one can show that 
two solutions of (8) are singular. Regularity at the horizon demands that these singular 
solutions should be absent. To specify Cauchy data at the horizon we must therefore 
specify just one more quantity at the horizon — conveniently ai(l). Also k must be 
specified. At infinity there are two exponentially growing solutions to (8). The natural 
boundary conditions at infinity are that these solutions must be absent. Thus we have 
two free parameters, ai(l) and k, and two boundary conditions at infinity to satisfy. An 
appropriate two-parameter shooting algorithm is easy to implement. The only subtlety 
worth mentioning is that because we have two growing solutions at infinity, it's possible 
to them to compete so that ai(y)^ + Ci(y)^ is small at a single large value of y, but 
large at other vahics. So to test whether one has eliminated both growing solutions, 
an efficient and robust method is to integrate ai{y)^ + Ci(|/)^ over some interval near 
the end of the integration region. We have done this, with the result ai(l) = —0.552 
and k — 0.876. This value of k is in accord with the results of [1].^ 

3.3 Mass, entropy, and temperature 

Before going on to O(A^), let us describe in more detail how the mass and entropy can 
be calculated. The results in this section are vahd at any order in A, and indeed make 
almost no reference to the expansion (5). Following [6], we define the mass as 



where is the trace of the extrinsic curvature of a surface E at constant y and t, 
computed with respect to the four-dimensional spatial metric: that is. 



^In comparing with Figure 1 of [1], one must note that D = 4 in [1] corresponds to the five- 
dimensional black string in our language. Also, it appears that the numerical integrations in [1] were 
performed at Schwarzschild radius r+ = 2 (see also the remarks on p. 17 of [2]). To convert to r+ = 1, 
and n should be doubled. And /i is what we have called k. Thus the intersection of the interpolating 
curve with the horizontal axis at « 0.43 in their Figure 1 is indeed in accord with k = 0.876. 



2ey^a, + (1 -Ay) a[ - 2 {c[ + (-1 +y) y<) = 
2ai+(2 + A;^ (3-4y) y^) ci + 3a[-3ya[ 

+9c[ - 16yc[ + 8y^ c[ - y{-l + y){3 - Ay) d[ = 



(8) 




(9) 




(10) 
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where ^gr is the metric on a. t — const spatial shce, and is the unit normal to 
the surface y — const in this slice. The quantity ^Kq is the extrinsic curvature for a 
reference geometry, in our case flat space. A large y limit is taken after subtracting the 
contribution of the reference geometry, and then the energy should be finite. In making 
the subtraction, the surface S must be chosen to have the same induced geometry in 
the black string background and the reference geometry, at least up to a sufficiently 
high order in large y. In our case, we may altogether ignore the functions for i 7^ 
(effectively setting them equal to 0) because at large y they fall off exponentially in 
y, whereas the terms in — ^Kq that make finite contributions to the energy are 
0(l/y^). Let us assume the leading order behavior 

B = ^ + 0{l/y'') C^^^^^ + 0{l/y). (11) 

y y 

Then, suppressing all exponentially small corrections, we can write 

= - - 2Coo^ - L(i + _ 2Coo) + 0(1/2/') , 

y y^ y^ 

This is to be compared with ^Kq — 2/yQ for the flat reference spacetime. However, in 
making the subtraction one does not want y^ = y, but rather yo = ye^ = y + C oo ^ogy + 
0{l/y) so that the S'^ has the same curvature in the black string geometry and the 
reference geometry. Thus we obtain 



y y (13) 

'X - = L(i + _ 2Coo) + 0(l/y') . 
y 

The area of the S'^ is 471?/^, and the length of the is 27r/A;, so we see that the energy 
of the black string is (1 + 3i?oo — 2Coo)/A; up to overall constants. In later sections we 
will allow k to vary, so it is more convenient to define e = kE/2n as the average energy 
density. Then 

- = 35^-20,0. (14) 

e 

In applying the result (14), it must be recalled that we assumed (11). In a perturbative 
treatment in A, one must check this asymptotic behavior at any given order. Most 
modifications of the asymptotics would render the mass infinite. 
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Computing the entropy is more straightforward: it is just the horizon area over 4G5, 
so the fractional change in the average entropy density, s — kS /27r, is 



Ss 



-1, (15) 



where () means to average over the x direction. Thus {5s)/s winds up being some 
combination of the Xi evaluated at y = 1. 

The temperature is also easy to compute, using for instance the standard prescription 
of rotating to Euclidean signature and demanding no conical deficit at the horizon. The 
result is 

^ = e^-^ - 1 (16) 

The right hand side is not manifestly independent of x. A convenient check on boundary 
conditions at the horizon is to verify that they preclude any contribution to (A — S)|^^^ 
other than from zero modes. 



3.4 The zero modes at 0{\^) 

As per the discussion of the previous section, to compute {Se)/e and {Ss)/s to O(A^), 
all we require is the functions Xq to this order (where as usual, X can be a, b, or c). 
So we postpone discussion of the X2 until section 3.5. 

The equations for the zero modes at O(A^) are rather complicated. All the relevant 
equations take the schematic form Xq = Xf, where the left hand side denotes some 
linear combination of ao, bo, cq and their first and second derivatives, and the right 
hand denotes a quadratic expression in ai, 61, ci and their first and second derivatives. 
Using (6), (7), and (8), one can eliminate all second derivatives from the right hand 
side of Xq — Xf, and also eliminate all dependence on bi and b[. (This is a good idea 
for numerics, because the functions Xi are known only numerically, and differentiating 
them leads to more numerical noise.) One of the resulting equations is 

-^2(-l + y)y-^bo = (^rf^(2 y {s - 7 y + Ay') + k' y (-3 + Ay) 

-2 i-l + y) ci {a[ + 2c[) 

-{3-5y + 2y') (2 (-1 + y) a[ - c[) {a[ + 2c[) 

+ ai (k'^y (-3 + 10y-8y2) ci + 2 (-1 + y) {a[ + 2c[))). 

(17) 

This equation is easy to integrate twice numerically. In order to have regularity at the 
horizon, &o(l) must be finite, and this fixes a constant in the first integration. The 
constant in the second integration can be fixed by demanding bo ^ as y ^ 00. This 
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boundary condition does not come from asymptotic flatness, but it could be arranged 
through some rcscalings. No shooting is necessary. After integrating, one can easily 
extract B^o ~ 0.41A^. 

The other equations derivable at 0{X^) for the zero modes take the form 

a'o - 4yao - 2co + 2yao -2y'^aQ = Si{ai,ci) 
-4co + 2a'Q - 2ya'Q + 6c'^ - 8yc', + 2yc'^ - 2y^ c'^ = S2{a,,c,;bo) (18) 
8 Co - S a'o + 8y aQ - 4cq + 8y Cq = S3{ai,ci]bo) , 

where the source terms 5*1, 5*2, and 5*3 are quadratic expressions in Oi, Ci and their 
first derivatives, like the right hand side of (17). 5*2 and also have /znear dependence 
on 6o- We regard 6o as fixed from the considerations of the previous paragraph. Note 
that ao appears only through its derivatives. An additive shift of corresponds to a 
multiplicative redefinition of t. The equations in (18) are not independent: the first, 
for example, follows from the second and the third together. For numerics, it would 
seem most straightforward to keep the third equation and drop either the first or the 
second, since the third is first order. However, it turns out to be better to drop the 
third equation and integrate the first two numerically.^ In order to do this we need 
Cauchy data at the horizon. If we assume regularity of ao, Cq and their first and second 
derivatives, then by setting y = 1 in the first and third equations in (18) we obtain 

-Sa'o - 2c'o = Si(ai, Ci) , -4co(l) - 2c^,(l) = ^2(01, Ci; 60) , (19) 

where all quantities arc evaluated at y = 1. The third equation in (18) does not lead 
to any additional boundary conditions at the horizon. Two more boundary conditions 
are required to complete the Cauchy data: one of them is the trivial additive constant 
on oo, which can be adjusted a posteriori so that ao(?/) — > as y — > 00; and the other 
is the value of co(l). Any choice of co(l) leads to a solution that is regular everywhere, 
with asymptotics as in (11). The O(A^) contribution to Coo depends on Co(l). 

Clearly this leaves us with a puzzle: why is Coo, and hence the mass, indeterminate? 
Roughly, the answer is that, at the order to which we are working, we are still free 
to "superpose" an arbitrary O(A^) change in the mass of the black hole on top of the 
change that the non-uniformity induces. The value of co(l) is undetermined at O(A^), 
and will be fixed, if at all, by considerations at higher order. (We will revisit this 
issue in section 3.5.) The equations (18) with the Si set to correspond precisely to 
perturbations of a Schwarzschild black hole in four dimensions in a gauge where the 
radial part of the metric is exactly dy'^/ (l — l/y). With this choice of gauge, an increase 
of the mass of the black hole is expressed as a solution to the homogeneous equations 

^T. Wiseman pointed out to me that this alternative led to considerably more tractable numerics 
near the horizon. 
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derived from (18) with Co{y) ~ C^{\ogy)/y for large y. Any multiple of this solution to 
the homogeneous equations can be added to a solution to the inhomogeneous equations: 
that is the origin of the ambiguity in the mass. 

To proceed, let us set co(l) = 0. Numerical integration of (18) is now straightforward. 
Summarizing the results of the numerics so far: 

ai(l) = -0.55 6i(l) = -0.55 Ci(l) = 1 

ao(l) = 0.53 6o(l) = 0.77 co(l) = (20) 
= 0.41A^ Coo = 0.29A2 . 

Given these values, together with (14) and (15), one may evaluate, to O(A^), 
5e 

— = - 2Coo ^ 0.64A2 

(21) 



5s_ 
s 



bo + 2co + |6i + 6ici + cl 



A^ ^ 1.29A2 . 

y=l 



Adjusting co(l) changes {Ss)/s by exactly twice the amount that it does {Se)/e: this 
is because the entropy-mass relation for the four-dimensional Schwarzschild solution is 
5" oc M^. A quantity which is supposed to be independent of co(l) is 

OS OG 

Aes^ 2— = (JiA^ where ai ^ 0.002 . (22) 

s e 

The coefficient ui is remarkably small compared to the values in (20). This may lead us 
to suspect that the true value is zero, so that the non-zero result in (22) is pure round- 
off error. We will presently give an entirely analytical demonstration that ai — 0; 
however, let us for the moment examine how we would judge just from the numerics 
whether (22) is consistent with a zero result. One way way to test this hypothesis is to 
vary Co(l) and see how much ai changes. Since all the values in (20) are on the order 
of unity, it would seem natural for Cq{1) to have a similar magnitude. So the deviations 
in the result that come from assigning Co(l) the values ±1 rather than are something 
hke a standard deviation. This "standard deviation" turns out to be roughly 0.02, so 
our result is indeed consistent with ai — (in fact, one might say that we were lucky 
to get a value as small as the one in (22)). Another test is to upgrade the numerics in 
various ways and see what changes. The results quoted in (20)-(22) come from the last 
of five iterations of improvement, which added about one decimal place of accuracy. 
The first three iterations were about one standard deviation positive. 

The bottom line is that numerics arc consistent with the hypothesis that, at O(A^), 
a uniform black string and a non-uniform one with the same mass also have the same 
entropy. More precisely, the non-uniform black string's entropy exceeds that of the 
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uniform black string only by about one standard deviation, as defined heuristically 
in the previous paragraph. Now let us show analytically that Ags = o(A^), using 
the First Law of thermodynamics (which certainly should hold for the non-uniform 
solutions at hand, since as static solutions to the vacuum Einstein's equation they 
have well-defined event horizons). We will need to assume that, if the asymptotic size 
of the is held fixed, then the mass is specified unambiguously once one has specified 
A. That this is true will emerge from the discussion in section 3.5. Thus there is a 
one-parameter curve of non-uniform solutions in the rj-X plane, emanating from the 
point (?7c) 0) ~ (0.876,0). Suppose we slide incrementally along this curve, away from 
{ric.O). Then the First Law says dE = TcdS, where Tc is the temperature at (r^c, 0). 
At this point, the black hole is uniform, so we may use the standard four-dimensional 
Schwarzschild relation = Ec/{2Sc)- Rearranging the first law slightly, we obtain 
2{dE)/Ec = {dS)/Sc. Thus indeed Aes = 0, at least to some lowest order. Now we 
have to explain why this works to O(A^). The reason is that the temperature deviates 
from Tc only at O(A^), and then this small deviation is further suppressed (to O(A^), 
in fact) when one integrates the infinitesimal form of the First Law out from to A. 

The assumption that the asymptotic size of the is constant is needed in the above 
derivation because otherwise the First Law would have an additional term correspond- 
ing to variations in this quantity. The assumption that only a one-parameter family 
of non-uniform solutions exists with fixed asymptotic size of is needed for the same 
reason: otherwise there would effectively be another thermodynamic observable that 
would contribute to the First Law. 

Obviously, although we may feel relieved that our numerics through O(A^) passes one 
non-trivial analytical check, we should also feel disappointed that no really meaningful 
numbers came out of the analysis so far, beyond what is in the literature. To do better 
we must proceed at least to O(A^). There are some conceptual subtleties involved, 
which we will discuss in the next section. 

3.5 Higher orders in perturbation theory 

To organize higher order in perturbation theory, it will be useful to consider the ex- 
pansions 

oo 

A^Yl ^""^n COS nKx 

n=0 
oo 

B^Yl ^"^n COS nKx (23) 

n=0 
oo 

C = ^ A"C„ cos nKx , 

n=0 
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where for X — A, B, and C we have the additional expansions 

oo 

X^ = Y.\'PX„,p, (24) 

p=0 

and now the X„ p are independent of A. Also K itself may be taken to have an expansion 
in A, 

oo 

K = J2\^%. (25) 

q=0 

Translating back into the notation of (5), we have ko = k, Aq^q = 0, Aq^i = qq, Ai^ = ai, 
A2,o = 0,2, and the same for B and C. We do not see any reason for non-analytic terms 
in A to appear at any order in (23) or (24), or for odd powers of A to appear in (24) 
or (25). If either thing happened, then it would seem hkely that there are additional 
physical constants of integration to specify beyond rj and the ratio of minimum and 
maximum radius for the horizon.^ 

Given the results of previous sections, it should be plausible that the Einstein equa- 
tions = boil down to an infinite set of ODE's of the form 

(26) 

The left hand side represents a first or second linear differential operator operating 
on the Xn^p. Linear equations for Xn,p arise at 0(A"+^^).^ For n 7^ 0, the form of these 
differential operators is identical to the forms in equations (6), (7), (8), only with k 
replaced by nko. The right hand side of (26) represents a sum of terms, each expressible 
as a product of several Xn^^p^ and their derivatives, as well as various kg., subject to 
two sum rules: J2i ^iUi — n, where each Sj = ±1, and J2i{ni + 2pi) + J2j qj — n + 2p. 
The first of these rules comes from Fourier analysis, and the second comes from power 
counting in A. This means that one may solve the equations (26) order by order in A 
(that is, in order of increasing n + 2p) — modulo a difficulty to be discussed below. At 
a given order in A, the equations with different n are independent. There is a further 
useful property of the equations (26), easily demonstrated order- by-order: for n ^ 0, 
one can find algebraic expressions for 5„p and B'^^ in terms of Cn,p, and other 
Xn,p of lower order, and first derivatives of these quantities. For n — 0, this is not 
possible. 

^Provided we take Ci.o(l) = 1, it is still tnio that the A of this section coincides with the one in 
(2), to leading order in A itself. Working out a more precise translation is straightforward once the 
Xn,p and kp are known up to a given order. 

^Starting at 0{X^), one also gets linear equations for functions Xn^p which were determined at 
lower order. These equations arise multiplied by powers of kp and x: that is, they are artifacts of 
expanding cosKx and sinKx in a basis of functions with period 2-K/ko. It was certainly easy to 
check at 0{\^) that these "secular" equations were automatically satisfied given the 0(A) equations; 
probably this could be demonstrated at all orders, but we will not attempt complete rigor here. 



12 



A subtler point is how we determine all the constants of integration, including the 
kq. There is some freedom on how this is done, corresponding to two obvious consid- 
erations: first, we have not entirely fixed the diffeomorphism freedom with the ansatz 
(4) and (23); second, starting at second order in O(A^), one has the freedom to simul- 
taneously add something to the mass and change K, in essence rescaling the solution 
by a perturbatively small amount. Even after demanding that A{x,y), B{x,y), and 
C{x,y) vanish in the limit y — > oo, there is some ambiguity left. To "fix a scheme," 
one could for example say Xq^q — 0, Co,p(l) = for all p, Ci,o(l) = 1, and Ci,p(l) = 
for p > 0. Then the kg would be fixed as an integration constant in the equations for 
Xi^q+i. We will call this the standard scheme. 

To make the discussion more definite, let us consider the standard scheme through 
the first few orders in perturbation theory. At 0(A), we have only the Xi^o equations, 
and these fix two integration constants, which in the standard scheme are ko and 
Aifi{l). At O(A^), the requirement that A{x,y) and C{x,y) should vanish as y — > oo 
fixes ^o,i(l) and So,i(l), while Co,i(l) is set to zero as part of the standard scheme. 
This is an arbitrary resolution of the ambiguity we encountered in section 3.4. Also 
wc have at O(A^) the equations for X2fi, and asymptotic flatness flxes the integration 
constants ^2,0(1) and C2,o(l)- Finally, at O(A^), we have the equations for Xi^i and 
X3 0- The latter are similar to the X2fi equations, and asymptotic flatness flxes A3 o(l) 
and C3 0(1); the former are similar to the Xi^ equations, and in our scheme of setting 
Ci,i(l) = they flx ki and Ai^i(l). As an illustration of how schemes might be altered, 
we could leave Co,i(l) free in the analysis of the Xq^i equations, but set ki — as well 
as Ci^i(l) = in the analysis of the 1 equations, which would leave us with Co,i(l) 
and 74i^i(l) as the constants of integration which are fixed by asymptotic flatness in 
the analysis of the Xi 1 equations. This scheme is computationally disadvantageous 
because when one changes Co,i(l), it is necessary to re-integrate the Xq^i equations 
before solving the Xi^i equations. Physically, the new scheme corresponds to holding 
fixed the asymptotic size of the around which the black string is wrapped, whereas 
the standard scheme allows this asymptotic size to vary at O(A^) but sets to zero 
a particular O(A^) contribution to the change in mass. Naively, one might propose 
yet another scheme where the total mass is held fixed at O(A^) and fci = as well, 
but Ci^i(l) is left free, in which case i(l) and Ci^i(l) would be the constants of 
integration fixed by analysis of the Xi.i equations. This seems natural, but in fact we 
suspect that this scheme is degenerate, because it takes us out toward non-zero A on 
a curve where r] is constant at least to O(A^). The dependence of 77 on A at lowest 
non-trivial order is part of the physical "output" of our numerics, and so should not 
be fixed as part of an arbitrary scheme. 

Let us now to summarize our immediate aims without getting too tangled up in 
schemes: we want to pin down the O(A^) change in the mass of the string which was 
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left undetermined in section 3.4. To do it we have to proceed to O(A^). A clever trick 
is to fix (5e)/e arbitrarily at O(A^), but instead allow the size of the to change at 
O(A^). That's equivalent because we can always rigidly rescale the whole solution to 
bring the size of the back to what it was, and in the process recover unambiguously 
the desired O(A^) contribution to 6e. 

So far, we have not discussed the X2,o equations in any detail; indeed, the analysis 
of them is a rather uninteresting replay of the Xi o equations. Let us only remark 
that B2fi can be determined algebraically, and the remaining two equations are second 
order in ^2,0 and 6*2,0 , and involve source terms which can be written as quadratic 
expressions in Ai.o, Ci,o, and their first derivatives. Regularity at the horizon fixes 
^20(1) and C2q{1) once ^2,0(1) and 6*2,0(1) are known. These latter two quantities 
must be fixed arbitrarily to carry out a numerical integration, and their values then are 
determined by the requirement that both /l2,o and 6*2,0 shrink exponentially for large 
y rather than growing. The result is 

A,o(l) = 0.34, 62,0(1) = -0.69. (27) 

The X2,o equations arc completely independent of the Xq^i equations, so the arbitrari- 
ness of 60,1(1) does not affect (27). 

At O(A^), the Xi,i equations and the ^3,0 equations are independent. We have not 
analyzed the X3,o equations in any detail, but clearly they will turn out just like the ^2,0 
equations: -83,0 will be determined algebraically, and As^i^l) and 63,0(1) will be fixed 
by normalizability at infinity. It is the equations that will primarily interest us for 
the rest of this section. Once again, Bi i can be algebraically eliminated, and what is 
left is linear second order equations for and 61,1, sourced by expressions up to cubic 
in the Xi,o, and also depending on ki, Xq^i, and ^2,0. The form of the homogeneous 
equations (that is, with the source terms set artificially to zero) is identicalto (18), only 
with ai replaced by and Ci by 61,1 (this is a general fact, indicated in (26) by the 
dependence of only on X and n, not p). In our standard scheme, ki and y4i,i(l) are 
determined by normalizability of Ai^i and 61,1 at infinity, so there is a two-parameter 
shooting problem. The results of numerics are 

fci = 0.70, A,i(l) = -0.24. (28) 

The best way to extract the effect on the mass is to return to our definition of rj. 
For the five-dimensional black string, the "average" Schwarzschild radius appearing in 
(1) is proportional to the average energy density e, so rj can be expressed as eK up 
to a factor which includes a power of Newton's constant. Combining the O(A^) effects 
from (21) and (28) (both of which were computed in the standard scheme), we find 

Sri Se SK , ^ 

— = — + = 771 A^ where 771 1.45 . (29) 

T] e K 



14 



Because rj is independent of rigid rescalings of the entire solution, rji should be scheme- 
independent. In fact, there is approximately 0.6% variation in r)i when Co,i(l) is 
changed from to 1, and much smaller variation when Ci,i(l) is changed from to 1. 

It is possible to extend the First Law argument of section 3.4 to a computation 
of the entropy difference of the uniform and non-uniform solutions to 0{\'^). Again 
holding the asymptotic size of the fixed, let us integrate the infinitesimal form of 
the First Law along the curve of non-uniform black brane solutions, starting from the 
point where it joins onto the uniform solutions, and assuming expansions 

oo oo oo 

M = J2MpX''P S = J2SpX''P T = Y.T,\^P. (30) 

p=0 p=0 p=0 

Starting with Mq = 2X05*0 (true because at A = the equation of state is the same as for 
the four-dimensional Schwarzschild black hole) one quickly obtains Mi/Mq = Si/{2Sq) 
and M2/M0 = S2/{2So) + Ti 5*1/ (4X06*0). Now let us inquire how the entropy of a 
uniform string changes if we change the mass by AM. Since S oc M^, we have 
(AS) /So = 2(AM)/Mo + [(AM)/Mo]^ What we really want is the difference between 
the entropy of a non-uniform string and a uniform one of the same mass. So we set 
AM = MiA^ + MaA^ and obtain, to O(A^), 

'S'non— uniform '^'uniform 'S'non— uniform '^'uniform 




So 

2AM /AM\ 



2 



TiMi /Ml 



(31) 



The first equality holds, even though ^uniform deviates from Sq by O(A^), because the 
numerator is an 0(A'^) quantity. In the final equality, we have used the fact that 77 
is proportional to M, and the constant of proportionality is fixed if the asymptotic 
size of S^ is held fixed. We have also used the facts that \^Ti/Tq = [ST)/T when the 
asymptotic size of is held fixed, and that {6T) /T — {SK) / K is a scheme-independent 
quantity at O(A^). The final two expressions in (31) involve only quantities computed 
at O(A^), so with (28) in hand one can compute, to the relevant orders, 

5T 5K ^^^^2 
-- — ^ -1.04A^ 

_ o (32) 

*Jnon— uniform '-'uniform \A ^ r> rr> 

= (72A where (T2 ~ —0.59 . 



'S'uniform 

The coefficients quoted were computed in the standard scheme. There is approximately 
3% variation in a2 when Co,i(l) is changed from to 1, and much smaller variation 
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when Ci^i(l) is changed from to 1. The value of (T2 could be used as a consistency 
check on numerics at O(A^), just as the near-vanishing of ai was used in section 3.4 as 
a consistency check at O(A^). 

The 0.6% variation in rji upon changing Co,i(l) from to 1, and the 3% variation 
in (T2, are indications of the size of numerical errors, since in principle the quantities 
in (32) are scheme- independent. Additively, the variation in (T2 is only twice as large 
as the variation in 771, and it is nearly the same as the variation in ai discussed in 
the paragraph following equation (22). The upshot is that the variations observed in 
these three quantities are consistent with one another, and small enough that we can 
be quite confident that 771 > and (T2 < 0. 

3.6 Improving and extending the numerics 

There are various ways in which the numerics discussed in previous sections could be 
improved and extended. To begin with a rather technical point, we have run all our 
numerics with Mathematica's built-in NDSolve routine, which uses adaptive step-sizing: 
the discrete steps can change each time one solves a differential equation with different 
initial conditions. This results in decreased stability for shooting algorithms which 
became particularly noticeable in the analysis of the X2.0 equations. Fixing the step 
size once and for all would have its advantages where stability is concerned, although 
one would probably want to retain the feature that shows up with NDSolve, that the 
steps are smaller near the horizon. Another technical point which could clearly stand 
some improvement is the treatment of initial conditions. The rough-and-ready method 
we have used is to observe that second derivative terms cancel out of the differential 
right at the horizon, so evaluating the equations there results in a constraint on first 
derivatives — but one which we then impose at a slight distance away from the horizon 
(on the order of 10""^) so that NDSolve won't be faced with excessively small coefficients 
on the second derivative terms. A better method, though more labor-intensive, would 
be to obtain analytic approximations to the regular solutions to some moderate order, 
and use them as the seed for numerics. (This is straightforward in principle because 
the horizon is a regular singular point of the differential equations). It might even be 
hoped that if one knew the near-horizon asymptotics as well as the long-distance tails 
to sufficiently high order, a reasonably uniform approximation to the Xn,p could be 
obtained via matching. This is not an approximation that could be controlled by a 
small parameter, as it is in black hole absorption at low energy. It might nevertheless 
be useful for large n, where numerics gets harder and harder due to the increasingly 
powerful exponentials in the growing solutions. 

To determine Coo at any even order in A, we have to distinguish a (logy)/y behavior 
from a 1/y behavior in a function, Co{y), known only numerically. This is one of the 
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hardest aspects of the numerics, so it is natural to suspect that it is a dominant source 
of numerical error. The quantities ai, rji, and (J2 all depend on evaluations of Coo, and 
all three displayed variations of similar magnitude (of order 10~^) when the scheme 
was changed. On the other hand, spot checks of the values of the components of the 
Ricci tensor gave values on the order of 10~^^ up to O(A^), when numerical evaluation 
was performed after elimination of all second derivatives. This is only ten times the 
working precision for most executions of NDSolve. These numerical observations are 
consistent with the hypothesis that the determination of Coo dominates the numerical 
error. ^'^ 

In the end, customized code would be needed to optimize results. We have aimed in 
this paper to avoid excessively hard-core numerics, extracting instead what information 
we could from standard tools. The total amount of CPU time needed to check all the 
numbers presented in this paper is roughly an hour on a 800 MHz Pentium PC. 

3.7 Interpreting the numerics 

Broadly speaking, the goal of this paper has been to make some dent in the prob- 
lem of constructing stationary non-uniform black brane solutions and elucidating their 
thermodynamic properties. Our point of entry was the observation that there is a 
zero-mode fluctuation for the uncharged black string in five dimensions, with wave- 
number k — 0.876/i?schwarzschiid- This suggests that there is a one-parameter family 
of non-uniform solutions in the 77-A plane that joins onto the uniform solutions at a 
"critical point," {rjc, Ac) = (0.876, 0). We may then ask, 

1. What curve in the r]-X plane is traced out by these non- uniform solutions? 

2. What is the entropy along this curve? 

3. How does a black string behave in real time as one slowly adds mass or allows 
mass to slowly Hawking radiate away, such that rj passes through r/c? 

A perturbative expansion in small A seems well-suited to study the first two questions 
near the critical point. Given the answers to these two questions, we can apply the 
Second Law plus some "natural" assumptions to give a qualitative answer to the third. 

^°In the original hep-th version of this paper, a different source of numerical error dominated: the 
numerics on the second and third equations of (18) was unstable near the horizon, leading to a poor 
value for ao(l). This is an instance where matching an analytic near- horizon solution to a numerical 
one would have avoided any problem. 
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Figure 1: (a) Possible curves, drawn schematically in the rj-X plane, describing station- 
ary, uncharged, non-uniform black strings in five dimensions. The solid curve is the 
likeliest, (b) a versus A for these non-uniform strings. Again, the solid curve is the 
likeliest. 

Possible answers we could envisage to questions 1 and 2 are sketched in figure 1. 
Given the analyticity properties in A of the perturbation expansion, we expect 



V = Vc + J2 ^^^Vp 
p=i 



S- 



(33) 



a 



noil— uniform 



uniform 



1 + E A^v^ 

p=l 



p ' 



for some coefficients r]p and ap. In the first line of (33), as in (31) and (32), 5'non-umform 
and S'uniform are the entropies of black string solutions with the same mass}^ The 
First Law guarantees that ai = 0, and this was verified to good accuracy by the 
numerics through O(A^) (see (22) and the discussion below it). Numerics through 
0{X^) led to a positive result for rji (see (29)). Combining these numerics with First 
Law considerations further predicts (72 < (see (32)). Thus by our calculations of the 
previous few sections, we rule out the curves marked (3) in figure 1. 

Although we have no direct evidence, it would seem sensible that rj < rjc and cr > 
1 for sufficiently large A. The first inequality guarantees that, for some range of rj 
less than rjc, there are stationary solutions with finite A and a single maximum (n = 
1). The second inequality says that these solutions are accessible as the endpoint of 
Gregory-Lafiamme evolution for an unstable uniform string of nearly the same mass. 



Again we are glossing over a difference in the simplest conceptual definition of A, namely (2), and 
the computationally convenient definition of A, which in the standard scheme could be expressed as 
the Fourier coefficient of cos kx in the fractional deviation of the local Schwarzschild radius from its 
mean value. We do not think this ambiguity will affect the issue of analyticity, but clearly it would 
change the values of the coefficients in (33). 
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Taking these inequalities as working assumptions, we rule out the curves marked (1) 
in figure 1(a) and 1(b), and settle on the curves marked (2) as the likeliest. We have 
incorporated an additional prejudice into curve (2) of figure 1(b): namely, we expect 
that the non-uniform black string with maximal rj should have an entropy less than a 
uniform black string of the same mass, so that if one started with such a non-uniform 
solution and added mass, it would have someplace to go (namely, to a uniform solution). 
We emphasize that our current numerics determines only the leading non-zero behavior 
of T) and a for small A. The rest of our expectations are motivated by the hope that 
once wc know the continuous moduli space of static non-uniform solutions, wc will be 
able to give a complete qualitative account of the real-time dynamics in the vicinity of 

If indeed the curves marked (2) are correct, then the transition between uniform 
and non-uniform strings is first order. (The transition is also first order if there is 
a discontinuous moduli space, as discussed in the previous footnote. We will not 
consider this option further since we have no way of exploring it). Suppose we start 
with a very massive uniform black string wrapped on an whose size we hold fixed. 
Then as the string Hawking radiates, r) decreases from its initial large value. When 
it reaches r]A, there is another solution available, but provided fluctuations are small 
and tunneling is suppressed {i.e. -RschwarzschiW ^ ^Pi) the system will be unable to 
jump from (A) to (B) even if (B) has higher entropy (which it well may not, since 
(72 < 0). Instead it will proceed down to (C), become locally unstable, and evolve in 
some non-adiabatic fashion, settling down presumably to (D), which is possible since 
we have assumed in figure 1(b) that (D) has larger entropy than (C).^^ Going the 
other way, if one starts with a non-uniform string with f] < r]c and slowly increases 
its mass (for instance by feeding it dust), then the system would reach point (B) and 
then evolve non-adiabatically to (A). Again this is possible since wc have assumed in 
figure 1(b) that the entropy at (B) is less than at (A). We are tempted to say that the 
latent heat of the transition from uniform to non-uniform solutions and back would 
be determined by a at points (B) and (D), but this seems a slight misnomer since in 
the usual language of phase transitions, both jumps (from (C) to (D) and from (B) to 
(A)) would be best described as discontinuous transitions from a meta-stable phase to 
a stable one. 

To test the picture suggested in the previous paragraph, it would be very interesting 
to go to O(A^) in the numerics and obtain the constants 772 and a^. A certain range 

^^Of course, one might also imagine the more arcane possibihty that curve (1) is right and that 
there is a new disconnected branch of non- uniform solutions at finite A which extends below r]c. 

^^Since there is some non-adiabatic behavior, a finite amount of mass will be lost to classical 
gravitational radiation. However, this lost mass will scale as some positive power of the characteristic 
frequency of the transition. By scaling up the size of the whole system, we can make this loss 
controUably small. 
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of negative values of r]2 and positive values of would support the hope that a < 1 
at (B) at cr > 1 at (D). Unfortunately, because (B) and (D) arc both at finite A, 
there is a possibility that low order perturbative results will mislead. Also, one would 
have to be careful to use an unambiguous definition in A. In computing lowest order 
results, it's OK that we used interchangeably two different definitions of A which differ 
multiphcatively by 1 + O(A^). That won't work at higher orders. 

One might imagine situations in which the curves marked (3) in figure 1(b) are 
qualitatively correct — for instance, in a different dimension, or with some form of 
charge on the string. Then the transition is probably continuous: a uniform black 
string crossing through r)c will smoothly develop non-uniformity by moving onto the 
branch of moduli space where A 7^ 0. The Second Law permits this provided a > 1, 
as we assumed in curve (3) of figure 1(b). If a continuous transition were found for 
some black brane system, it would be interesting to compute critical exponents, though 
naively it seems pretty likely that everything of interest will be analytic in A, perhaps 
even A^. 

Another way to describe the difference between a first order and continuous tran- 
sition, as we are using the terms, is that in a continuous transition, the amount of 
"thrashing," or non-adiabatic evolution, that an initially uniform string below rjc expe- 
riences in its real time evolution before settling down to a static non-uniform solution, 
would vanish as the initial rj approaches r]c] whereas in a first order transition there is 
a finite minimum amount of thrashing. 

Stranger possibilities than curves (2) or (3) may be contemplated: for instance, if 
the curves marked (1) are correct, then the non- uniform solutions in the moduli space 
we have explored are never reached starting from a Gregory- Laflamme instability, both 
because they have the wrong mass and because their entropy is too low. This scenario 
is not ruled out by our numerical results, but we would regard it as very surprising. It's 
worth noting that real-time simulations of gravitational collapse could lead to static 
non-uniform solutions regardless of which curves arc right in figure 1. 

Clearly we are in the realm of the speculative, and would be considerably aided by 
results at higher order in A. Some form of real time numerics would also be desirable in 
order to get to finite A. An interesting point which might be cleared up by analytical 
means is whether one can show that non- uniform strings do not exist above a certain r). 
It seems likely that this is true, since for very large 77 we could perform a Kaluza-Klein 
reduction along the and then try to argue that static spherically symmetric black 
holes in the resulting four- dimensional theory are unique. 

^*We thank G. Horowitz for raising this point. Note, this hne of argument might be more subtle 
than it naively appears, since if static non-uniform solutions exist at all in five dimensions, the four- 
dimensional uniqueness argument cannot be true for all 77, but must somehow hinge on the relative 
size of typical Kaluza-Klein masses and the inverse Schwarzschild radius. 
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4 Other black branes 



Although relatively unexplored, the phenomena associated with non-uniform black 
branes promise to be quite interesting and varied. The main difficulty is the dearth of 
problems that can be attacked without the help of heavy-duty numerics. Thus, at this 
stage, we have many more questions than answers. 

The main question we have addressed in this paper is the nature of the transition 
from uniform to non-uniform solutions at the critical size where the Gregory-Laflamme 
instability first shows up. Focusing on the uncharged black string in five dimensions, 
we have argued that although there is a continuous moduli space of solutions, with a 
branch of non-uniform solutions merging onto the uniform ones at 1] = r^c, the shape 
of the curve in the rj-X plane is probably such that there is a ffist order transition. 

Since our calculation was clearly quite specific to black strings in five dimensions, let 
us now consider what our expectations might be for more general black branes, based on 
symmetries and genericity. Before we begin, it is worth considering the analogy with the 
Landau treatment. Generally speaking, crystallization is a first order transition, and 
the reason is that there are cubic terms in the Landau free energy. More explicitly (see 
for example [7] for even more detail) one imagines starting with a uniform phase of some 
substance which at some temperature manifests an instability in density perturbations 
at a particular wavelength, call it k. We might write the density perturbations in 
terms of a scalar field, ~ 6p. Then the instability is toward developing some VEV 

— * 

for (f) — (f)j^cos{k • x). We might suppose that the instability in develops in a smooth 
fashion — for instance, via terms in the Landau free energy of the form 



where -u > and r(T) passes through zero at the crystallization temperature. If (34) 
were all, then the transition would be second order. But translation invariance also 
allows cubic terms of the form 



Such terms are allowed provided there is no —> —0 symmetry, and provided that we 

— * — * — * 

can construct equilateral triangles of allowed momenta, ki, k2, and k^. If the cubic 
terms (35) are present, then inevitably there is a first order phase transition to a state 
where several different wave-numbers condense simultaneously, with phases set in such 
a way as to make ^3 < 0. We will refer to this as "phase locking." In fact, one might 

^^This analogy was developed in part through conversations with participants of "Avatars of M- 
theory" and subsequently with E. Witten. 



(34) 
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further expect that the preferred crystal structure in dimensions higher than two (where 
there is only one natural candidate emerging from (35), namely the hexagonal lattice) 
will be the one with the most equilateral triangles per unit cell. This actually tallies 
with reality in three dimensions for crystal structures near their melting point: a great 
many are BCC 

The non-uniform black hole problem is similar to the model of crystallization de- 
scribed above in two important respects: 

1. Non-zero wavelength modes are unstable. 

2. There is nothing hke a ^ —cj) symmetry, since it's obviously different to make 
the black hole horizon bigger than to make it smaller. 

But it's different in three important respects: 

3. All wavenumbers below a critical one are unstable. 

4. We must work in the microcanonical ensemble. 

5. For large, smooth, classical solutions, no tunnehng is allowed. 

Compactifying the dimensions parallel to the black hole horizon is a convenient way of 
controlling point 3): we can contrive to deal with only one unstable mode at first (just 
below 77c), then try to work up to several unstable modes at smaller 77. Compactifica- 
tion has the additional advantage that one ducks out of any Coleman-Mermin- Wagner 
worries about whether long-range order is possible in low dimensions. 

Based on the computations in section 3, it appears that the black string transition is 
first order, but for reasons having nothing to do with cubic terms like (35), which are 
of course impossible in one spatial dimension. A closer analogy in a Landau treatment 
would be a free energy like = r(T)0^-|-M0*^+w0®, where u is negative, v is positive, and 
r(T) varies. For such a free energy, if tunneling and bubble nucleation is suppressed, a 
system that starts at = for positive r{T) stays there until r{T) — 0, then rolls off 
to a lower minimum at finite 0. For black strings in higher dimensions, or for charged 
black strings, we would not be at all surprised to find second order transitions. The 
positivity of rji, like the negativity of u in the analogy above, seems more likely to be 
a fact specific to the five- dimensional black string than a general property of all black 
brancs experiencing a Gregory-Lafiamme instability. 

For p-branes with p > 1 , can we infer that transitions from uniform to non-uniform 
branes should be generically first order? The answer depends on what one really means 
by a first order transition for black branes. What does seem likely is that the moduh 
space of non-uniform branes for p > 1 will extend both above and below 77 = 77c, due to 
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phase-locking effects. One could then imagine tunneling from uniform to non-uniform 

branes at some r) > rjc — that is, before the uniform brancs become unstable. This would 
be the closest analogy to what is usually meant by a first order transition in condensed 
matter physics. However, if we restrict ourselves to systems which evolve classically 
and have only small fluctuations, then tunneling is suppressed, and one can only ask 
what happens to a uniform brane once it starts experiencing the Gregory-Laflamme 
instability. As discussed above, our preferred notion of a first order transition in this 
setting is a transition that takes place non-adiabatically and results in a finite minimum 
non-uniformity. Then what determines whether a transition is first order or continuous 
is whether rj > rjc or r] < rjc for small non-zero A. In a perturbativc treatment, we expect 
that the leading contribution to rj will generically be at O(A^). For p > 1, it seems 
likely that there are odd powers of A contributing to rj; however we cannot see how an 
0{X) contribution would arise: such contributions, it seems, could be soaked up into 
a redefinition of rjc- So a determination of the O(A^) contribution to rj, similar to the 
computation of rji in (29), is the crucial indicator of first or second order behavior. If 
this contribution is positive, then a curve like (2) of figure 1 would pertain, and the 
transition is first order; if negative, then a curve like (3) would pertain, at least for 
small A, and the transition is continuous. 

Of course, if one is in a regime where tunneling is appreciable, then what's most 
relevant in all cases (including the black string) is whether there are, at a given value 
of 77, more entropic solutions at any finite A than the uniform black string. But in such 
a regime, other effects, like Hawking radiation and the classical gravitational radiation 
during non-adiabatic evolution, also become appreciable. 

This paper has been focused largely on 77 close to rjc — that is, on black strings com- 
pactified on circles on the order of their Schwarzschild radius. Let us now turn our 
attention to small 77, i.e. the limit where the black string or brane is decompactified. 
Then there would be many unstable modes on a uniform black string or brane (in- 
finitely many in the strict 77 — > limit). But, as is evident from Figure 1 of [1], the 
most unstable mode (that is, the one with maximum Q) is at approximately k — kc/2, 
where kc corresponds to the shortest wavelength that is unstable. The dispersion re- 
lation for the unstable modes, as computed numerically in [1], is (very approximately) 
n fa ^sinnk/kc. Thus if one starts with a uniform black string with very small 77, 
the "dominant" instability has finite wavelength, and one might optimistically expect 
a crystal structure to form with lattice spacing a ^ An/kc. A serious issue, however, 
is whether any crystal structure would be stable. Intuitively speaking, there are two 
potential threats to stabihty: 

1. Long- wavelength phonons: at wavelengths much long than Air/kc, it might seem 
plausible that something like the Gregory-Laflamme instability still operates, 
tending to cluster more mass in one region (populated by many maxima) at 
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the expense of another. Actually, such an effect might even manifest itself at 
V ^ where we could start with a black string (say) with two identical 

maxima and see whether the system is stable toward one maximum growing 
while the other shrinks. 

2. Short- wavelength instabihties: if the "necks" between maxima get too long and 
slender, then the local scale of Gregory-Lafiamme instabilities might be pushed 
small enough to occur locally in the neck. This might for instance create new 
miniature local maxima in the nooks between the big ones. 

By focusing on solutions with 77 close to rjc, we have avoided the long-wavelength 
phonon problem entirely: they are projected out by finite volume. As long as A is less 
than about 1, we can also be fairly confident that short-wave instabilities don't crop 
up. However, if there really is first order behavior in the five- dimensional black string, 
rather than a continuous transition at 77 = r]c, then we cannot completely rule out the 
possibility that short-wavelength effects will destroy the stability of the n = 1 solutions. 
The criterion for these UV instabilities to show up is a region of the black brane which 
is considerably less "thick" (as measured by its average Schwarzschild radius) than it 
is wide (as measured in the directions parallel to the horizon). 

If the instabilities described above are absent and there is a genuinely stable crystal 
structure at 77 = 0, then of course the physics at finite rj would be largely determined 
by frustration: if a black brane is compactificd on a torus which is commensurate with 
its stable crystal structure, then the torus is just filled with a region of undcformcd 
crystal; whereas for other sizes or shapes, the system fits in about as many unit cells as 
it can, with some strain because of the boundary conditions, and generically first order 
transitions between n maxima and n + 1. The simplest reason that such transitions 
should be first order is that they have different symmetry groups, neither of which is a 
subgroup of the other. For instance, in the case of a black string, the symmetry group 
of a solution with n identical maxima is a semi-product of the cyclic symmetry group 
Z„ and the reficction symmetry group Z2. 

At this point, we do not see a truly compelling reason to believe that the instabilities 
described above are absent at small rj for the generic non-uniform black brane solution. 
If they are present, then there would be no foreseeable endpoint to the evolution of an 
unstable black brane that extends over a volume much greater than its Schwarzschild 
radius. The evolution might, for instance, pass "fairly close" to a crystal structure, 
but then be driven away by a long-wavelength phonon instability. It would be difficult 
to check whether this happens via straightforward numerical solution of Einstein's 
equations in real time. The reason is that one needs an arbitrarily large range of the 
coordinate we have called x in order to include the effects of soft phonons. If numerics 
show that stable solutions form at 77 = 0.977c, and at 77 = 0.577c, and at 77 = O.I77C, a 
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skeptic might still point out that ai rj — O.Olrjc there are ten times as many phonon 
modes available, and the softest could be tachyonic. Perhaps one could work at a 
series of finite but small values of f] and then try to extrapolate the observed phonon 
dispersion relation uj{k) down to zero wavenumber. The short-wavelength instabilities 
seem somewhat less of a worry: if they were shown to be absent for finite r] (say in a 
real-time numerical treatment of n = 1 solutions) then it would be a surprise to find 
them at small 77 unless triggered somehow by a soft phonon instability. 

To improve the analogy with crystallization, one could consider a Landau free en- 
ergy in which all sufficiently soft modes of are unstable at quadratic order. What 
structures form and are stable is then highly dependent on the structure of the 0^ 
and 0^ terms. One could contrive, for instance, for a VEV of the most quadratically 
unstable mode to stabilize all the other modes — or for such a VEV to actually destabi- 
lize modes which were stable at = 0. The possibilities for non- uniform black branes 
seem, at this point, potentially as varied and interesting. It would be a fascinating 
enterprise (though, most likely, a computationally intensive one) to investigate which 
of the possibilities are realized by simple gravitational lagrangians. 

Some final remarks are in order regarding the case of charged black branes. For a 
black p-brane, we wish to consider only charge under a. p + 1-form gauge potential. 
The charge determines a definite length scale -Rchargc to which which the horizon radius 
-^horizon should bc Compared. If i?horizon ^ -Rchargc, then the physics should be essentially 
the same as for the uncharged case. When i?horizon ~ -^charge; the black brane often 
becomes locally thermodynamically stable. It is has been argued in [8, 9] that for 
locally thermodynamically stable branes, there is no perturbative Gregory- Lafiamme 
instability, no matter how small 77 may be; and it was further conjectured that the 
Gregory-Laflamme instability arises precisely when local thermodynamic stability is 
lost. Some evidence for this claim was presented in [8, 9], and further arguments have 
appeared in [10, 11]. A close look at these papers suggests that when the instability 
first arises (that is, when the black branes are just massive enough to be unstable), the 
wavelength of the unstable modes is extremely large. The short- wavelength instabihties 
described above are less of a worry: regions where the mass is decreased will eventually 
fall into thermodynamic stability, and then it can't have any local instabilities at all. 

It is perhaps worth reviewing the argument [8, 9] for why local thermodynamic sta- 
bility precludes a perturbative Gregory- Lafiamme instability. We emphasize "local" 
and "perturbative" here because it is always possible to increase the entropy of an 
infinitely extended, near-extremal charged brane by moving all the non-extremal mass 
into enormous, sparsely distributed Schwarzschild black holes through which the other- 
wise extremal charged brane now passes. The essence of local thermodynamic stability 
(which for a p- brane with no other quantum numbers than its mass and charge is sim- 
ply the positivity of the specific heat) is that if one keeps average mass density constant 
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but tries to make the system slightly non-uniform, the entropy goes down. Thus from 

the often-invoked Second Law, we learn that a Grcgory-Laflamme instability is not 
possible. The point of the conjecture [8, 9] that a Grcgory-Laflamme instability arises 
as soon as local thermodynamic stability is lost is that an infrared instability would 
seem plausibly to be entirely driven by thermodynamic considerations. 

5 Conclusions 

Although there seems to be a continuous moduli space of non-uniform black branes 
connecting to the uniform ones at the critical mass density where a Gregory- Laflamme 
density sets in, we have seen that this is not enough by itself to guarantee a continuous 
phase transition from uniform to non-uniform solutions. The results of our numerical 
study essentially boil down to two numbers, rji and a2- The positivity of rji (see 
(29)) means that slightly non-uniform solutions are above the critical mass density, by 
an amount proportional to the square of the amplitude of the horizon fluctuations, 
which we have called A. The negativity of a2 (see (32)) means that these non- uniform 
solutions have /oioer entropy than uniform black strings of the same mass, by an amount 
proportional to the fourth power of A. These results are exactly contrary to our original 
expectations, which were based on the hope that there would be a continuous transition 
describable as motion along the moduli space. The simplest scenario consistent with 
our numerical results, together with the general expectation [3] that the evolution of 
unstable branes settles down to a static endpoint, is that the moduli space of non- 
uniform solutions eventually curves down in the rj-X plane so as to provide the desired 
static endpoint, but at finite A; and that for the non-uniform solutions which are 
below the critical mass density, the entropy is higher than that of a uniform string 
of the same mass (see figure 1). If all this is true, then as soon as the mass of a 
five-dimensional, uncharged black string falls below a critical threshold, the Grcgory- 
Laflamme instability leads to non-adiabatic evolution to a finitely non-uniform solution. 
This is what we regard as a first order transition. 

How trustworthy arc the numerics? Clearly they are not perfect, since rii and o"2 
fiuctuate, respectively, by 0.6% and 3% as we "change the scheme" by an 0(1) quantity 
(as described more precisely in section 3.5). Abstractly, changing scheme means that 
we make some combination of a rescaling of the solution and a diffeomorphism on 
the X and y coordinates. The quantities rji and (72 are supposed to be scaling- and 
diffeomorphism-independent. It is quite plausible that the percent-level fiuctuations 
we observe in rji and (T2 are due to numerical error, particularly in the evaluation of 
the quantity we called Coo- We can say with some confidence that we got the signs 
right on rji and (J2. Since the claim of a first order rather than continuous transition 
for the five-dimensional black string is dependent only on the signs, we believe this is 
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a robust result. 

Extending the numerics to higher orders is a project that we hope to report on in 
the future. We also hope that it may become possible to address some of the issues of 
stability raised in section 4 at a sufficiently high order in perturbation theory. 
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